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Abstract. We discuss the structure of the equilibrium states of a regularized Keller-Segel model describing 
the chemotaxis of bacterial populations. We consider the limit of high degradation of the secreted chemical 
where analytical results can be obtained. Below a critical effective temperature, the system experiences 
a second order phase transition from a homogeneous phase to an inhomogeneous phase formed by two 
domains with uniform concentration separated by a thin interface (domain wall) . We study the properties 
of the interface and determine the bifurcation between a circular shape (spot) and a stripe as a function 
of the control parameters. We show the analogy with the structure of Jupiter's Great red spot which also 
consists of two phases with uniform potential vorticity separated by a thin annular jet. 
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1 Introduction 



The name chemotaxis refers to the motion of organisms in- 
duced by chemical signals p. In some cases, the biological 
organisms secrete a substance (pheromone, smell, food, ...) 
that has an attractive effect on the organisms themselves. 
Therefore, in addition to their diffusive motion, they move 
systematically along the gradient of concentration of the 
chemical they secrete (chemotactic flux). When attrac- 
tion prevails over diffusion, the chemotaxis can trigger a 
self-accelerating process until a point at which aggrega- 
tion takes place. This is the case for the slime mold Dic- 
tyostelium Discoideum and for the bacteria Escherichia 
coli. A model of slime mold aggregation has been intro- 
duced by Patlak [2] and Keller & Segel [3] in the form of 
two coupled differential equations. A simplified version of 
this model has been extensively studied in the case where 
the degradation of the secreted chemical can be neglected. 
In that case, the Keller-Segel equations become isomor- 
phic to the Smoluchowski-Poisson system describing self- 
gravitating Brownian particles |3j . An analytical study of 
this system of equations has been performed by Chavanis 
& Sire 5,6 7 8,9, 10 I11II12| in a series of papers . In partic- 
ular, self-similar or quasi self-similar solutions describing 
the chemotactic collapse have been obtained and the for- 
mation of a Dirac peak has been found in the post-collapse 
regime in d > 2 and in the collapse regime in d = 2. In 
parallel, a vast number of rigorous results concerning the 
existence and unicity of solutions of the Keller-Segel model 
and the conditions of blow-up have been obtained by ap- 



plied mathematicians (we refer to 13 for a connection to 
the mathematical literature). 

In this paper, we consider novel aspects of the Keller- 
Segel model. We first introduce a regularized model that 
prevents finite time-blow up and the formation of (unphys- 
ical) singularities like infinite density profiles and Dirac 
peaks. In this model, the local density of cells is bounded 
by a maximum value p(r) < a® which takes into account 
finite size effects and filling factors. Therefore, the Dirac 
peaks (singularities) are replaced by smooth density pro- 
files (clumps). With this regularization, there exists steady 
solutions (similar to the Fermi-Dirac distribution) for any 
value of the control parameter while the usual Keller-Segel 
model blows up above a critical mass M > M c (in dimen- 
sion d > 2). In addition, we consider the limit of high 
degradation of the chemical. In this limit, we show that, 
for sufficiently small (effective) temperatures T < T c , the 
system undergoes a second order phase transition from a 
homogeneous phase to an inhomogeneous phase. The bac- 
teria organize in two domains with uniform density p± 
separated by a thin interface. The resulting structure is 
similar to a "domain wall" in phase ordering kinetics |14j . 
We study in detail the structure of the interface (profile, 
width, surface tension,...) and determine the conditions 
for the bifurcation between a circular domain (spot) and 
a stripe in a square domain. This study can be performed 
analytically in the limit of high degradation (furthermore, 
our study is exact in d = 1). The case of a finite degrada- 
tion rate will be treated numerically in another paper. 

In previous papers | 15lHlll6 |. we have found a number 
of analogies between the chemotactic problem and other 
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systems of physical interest (self-gravitating systems, 2D 
vortices, Bose-Einstein condensation, Burgers equation). 
As mentioned previously, the Keller-Segel model is iso- 
morphic to the Smoluchowski-Poisson system describing 
self-gravitating Brownian particles jS] . In this analogy, the 
concentration of the chemical produced by the bacteria 
plays a role similar to the gravitational potential in as- 
trophysical systems (they are both solution of a Poisson 
equation) so that a number of analogies between biology 
and gravity can be developped 0]. In addition, the collapse 
of bacterial populations for M > M c or the collapse of self- 
gravitating Brownian particles for T < T c is by many re- 
spects similar to the Bose-Einstein condensation in phase 
space |T5]. Finally, there exists some analogies between the 
chemotactic aggregation of bacteria and the formation of 
large-scale vortices in 2D turbulence |15lH| . In that case, 
the concentration of bacteria plays the role of the vor- 
ticity and the concentration of the chemical produced by 
the bacteria plays the role of the streamfunction. In two- 
dimensional hydrodynamics, the vorticity field which is 
solution of the 2D Euler equation can achieve a statis- 
tical equilibrium state (on the coarse-grained scale) as a 
result of turbulent mixing (violent relaxation) |17j . In the 
two-levels approximation, the equilibrium vorticity pro- 
file is given by a Fermi-Dirac-like distribution [181119 ,20 
|2~T] . Interestingly, this is similar to the steady state of the 
regularized Keller-Segel model introduced in this paper. 
Furthermore, in the quasi-geostrophic (Q.G.) approxima- 
tion relevant to geophysical flows [22], the finite value of 
the Rossby deformation radius introduces a shielding of 
the interaction between vortices which is formally simi- 
lar to the degradation of the chemical in the chemotactic 
problem. In particular, the degradation rate k plays the 
same role as the inverse of the deformation radius R~ x . 
In the context of jovian vortices, Sommeria et al. \2'A\ and 
Bouchet & Sommeria [21] have considered the limit of a 
small deformation radius R — ► to account for the annular 
jet structure of Jupiter's great red spot. As we shall see, 
this is similar to considering a limit of high degradation in 
the chemotactic problem. Therefore, many interesting re- 
sults can be obtained by developing the analogies between 
these different topics. 

The paper is organized as follows. In Sec. El we in- 
troduce a regularized Keller-Segel model of chemotactic 
aggregation. We first provide a phenomenological deriva- 
tion of this model followed by a more kinetic approach. 
In Sec. |3] we study the equilibrium states of this model 
in a limit of high degradation. For T < T Cl we show that 
the solutions are formed by two phases in contact sepa- 
rated by an interface (the stability of the uniform phase is 
considered in Appendix EJ>. In Sees. 13. 113.31 we develop a 
"domain wall" theory to study the properties of the inter- 
face and determine its main characteristics (profile, width, 
surface tension,...). Asymptotic behaviors of these expres- 
sions are obtained for T — > T c and T — * in Sees. 13.41 and 
13.51 Analytical approximations of the jet profile are given 
in Sees. 13.61 and 13.71 in the form of self-similar solutions. 
Other approximations are given in Sec. 13.81 using match 
asymptotics. In Sec. 13.91 we show that the curvature ra- 



dius is constant so that, in two dimensions, the interface is 
either a line (stripe) or a circle (spot) . These results can be 
obtained equivalently by minimizing the free energy func- 
tional associated with the regularized Keller-Segel model 
(see Sec. I3.10jl . In Sec. 13. Ill we determine the phase dia- 
gram of the system and the range of control parameters 
(_B,T) where the equilibrium state is a stripe or a spot 
(the parameter B is related to the total mass of the con- 
figuration). Finally, in Sec. 0] we develop the close anal- 
ogy between our biological system and the jet structure of 
Jupiter's great red spot and other jovian vortices. 

2 The regularized Keller-Segel model 
2.1 The dynamical equations 

The general Keller-Segel model describing the chemo- 
taxis of bacterial populations consists in two coupled dif- 
ferential equations 

^ = V-( J D 2 Vp)-V-(£>iVc), (1) 



£^ = ~k(c)c + f(c)p + D c Ac, (2) 
at 

that govern the evolution of the density of bacteria p(r, t) 
and the evolution of the secreted chemical c(r, t) . The bac- 
teria diffuse with a diffusion coefficient D 2 and they also 
move in a direction of a positive gradient of the chemi- 
cal (chemotactic drift). The coefficient D\ is a measure of 
the strength of the influence of the chemical gradient on 
the flow of bacteria. On the other hand, the chemical is 
produced by the bacteria with a rate /(c) and is degraded 
with a rate fc(c). It also diffuses with a diffusion coefficient 
D c . In the general Keller-Segel model, D\ = Di(p, c) and 
D 2 = D 2 (p, c) can both depend on the concentration of 
the bacteria and of the chemical. This takes into account 
microscopic constraints, like close-packing effects, that can 
hinder the movement of bacteria. 

A very much studied version of the Keller-Segel model 
is provided by the system of equations 

|j = V • (DVp - xpVc) , (3) 
dc 

e-^ = D' Ac + ap- be, (4) 
at 

where the parameters are positive constants. Equation 
can be viewed as a mean-field Fokker-Planck equation as- 
sociated with a Langevin dynamics of the form 

civ , 

-= X Vc+V2DR(t), (5) 

where R(i) is a white noise and \ pl avs the role of a 
mobility. The stationary solution of Eq. Q is given by 

p = Aei c . (6) 
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This is similar to the Boltzmann distribution for a system 
in a potential — c. This suggests to introducing an effective 
temperature through the relation T e ff = D/x which is 
similar to the Einstein relation. For e = 0, the system 
©-© monotonically decreases (F < 0) the Lyapunov 
functional 



where C (p) is a convex function defined by 



1 f D 
F = — — / pcdr-\ / plnpdr, 

2 J X 



(7) 



which is similar to a free energy F = E—T e ffS where E = 
~ \ I P c( ^ v is the energy of interaction and S = — J p In pdr 
is the Boltzmann entropic functional 1 . For e ^ 0, the 
Lyapunov functional is 

F = ^- J [D'{\7c) 2 +bc 2 ] dr- J pcdv+— J plnpdr. 

(8) 

We shall consider here a more general situation where 
the mobility and the diffusion coefficient in the Keller- 
Segel model can depend on the density of bacteria. In 
that case, Eq. (0) is replaced by 



= V ■ MD(p)p) - X (p)pVc] 



(9) 



This can be viewed as a nonlinear mean-field Fokker-Planck 
equation |15j . It is associated with a Langevin equation of 
the form 



dr 

dt 



= X (p)Vc+^/2D(p)R(t). 



We define the functions h and g by 



dp 



xg(p) = px(p), 



(10) 



(ii) 



(12) 



where D and x are positive coefficients. With these nota- 
tions, Eq. © can be rewritten 



^=V.[Dh(p)Vp- X g(p)Vc}. 

Setting j3 = 1/T e ff = x/D, we obtain 
dp 



dt 



= V • [D (h(p)Vp - f3g(p)Vc 



(13) 



(14) 



This type of nonlinear mean-field Fokker-Planck equations 
has been discussed in They are associated with gen- 
eralized entropic functionals of the form 



5 : 



C(p) dr, 



(15) 



1 Note that these analogies with thermodynamics take even 
more sense if we remark that the Keller-Segel model is isomor- 
phic to the Smoluchowski-Poisson system for self-gravitating 
Brownian particles [1]. 



C"(p) 



Hp) 

9(p)' 



(16) 



The Keller-Segel model ©-® is known to exhibit 
blow-up solutions when the chemotactic attraction pre- 
vails over diffusion 13,10 . This reproduces the chemo- 
tactic aggregation of bacterial populations. In theory, the 
density can take arbitrarily large values and ultimately 
forms a Dirac peak. In reality, this singular evolution is un- 
physical as we expect finite size effects and close-packing 
effects to become important when the system aggregates 
and becomes dense enough. We shall regularize the prob- 
lem by introducing a sort of filling factor in the drift- 
diffusion equation {H}. Thus, we take h(p) = 1 and g(p) = 
p(l — p/<tq) so that Eq. (|3J is replaced by 



dp 
dt 



= V>[D (Vp-(3p(l-p/* )Vc)}. 



(17) 



With this modification, the mobility is reduced when the 
density becomes high enough (i.e. when p approaches the 
value (To) and this prevents singularities to form. Indeed, 
it can be shown that the density remains always bounded: 
p(r, t) < (To- This bound is similar to the Pauli exclu- 
sion principle in quantum mechanics, but it occurs here 
in physical space. The regularized drift-diffusion equation 
(|17fl was introduced phcnomcnologicallv in 15,4 to avoid 
infinite values of the density. 

In the following, we shall take e = for simplicity 
(although many results which are valid at equilibrium are 
independent on this assumption). This is valid in a limit 
of high diffusivity of the chemical 25 . If we introduce 
the notations k 2 = b/D' and A = a/D', we obtain the 
regularized Keller-Segel model 



dt 



[D (Vp-(3p(l-p/a )Vc 



(18) 



Ac - k 2 c = -Xp. (19) 
The stationary solution of Eq. i|18|) is given by 



P = 



1 + e -P c + a : 



(20) 



which is similar to the Fermi-Dirac distribution in physical 
space. From this expression, we clearly have p{r) < <7q. 
Furthermore, the Lyapunov functional can be written in 
the form of a free energy F = E — T e ffS where 

E=-±Jpcdr = -±J [(Vc) 2 + k 2 c 2 } dr, (21) 

is the energy of interaction and 



S=-a liJLlnJL + il-JL) ln ( 1 -JL ]>dr , 



Co c 



co. 



(22) 
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is the Fermi-Dirac entropic functional in physical space. 
The distribution H2()J) extremizes the free energy at fixed 
mass. Indeed, writing the first order variations in the form 
SF + aT e ffSM = where a is a Lagrange multiplier, we 
recover Eq. (|2U|) . Furthermore, it can be shown that a 
stationary solution of Eqs. (|18|l - l|19|l is linearly stable if, 
and only if, it is a minimum of F at fixed mass |15| . 



2.2 Phenomenological derivation of the model 

In this section, we develop a connection between the chemo- 
tactic problem and thermodynamics. To that purpose, we 
introduce the entropic functional (|22|l from a combinato- 
rial analysis which respects an exclusion principle in phys- 
ical space. Then, we obtain the dynamical equation l|18(l 
from arguments similar to the linear thermodynamics of 
Onsager. 

We divide the domain into a very large number of mi- 
crocells with size h. We assume that the size h is of the 
order of the size of a particle so that a microcell is oc- 
cupied either by or 1 particle. This is how the exclu- 
sion principle is introduced in the problem. We shall now 
group these microcells into macrocells each of which con- 
tains many microcells but remains nevertheless small com- 
pared to the spatial extension of the whole system. We call 
v the number of microcells in a macrocell. Consider the 
configuration {rii} where there are n\ particles in the 1 st 
macrocell, in the 2 nd macrocell etc., each occupying 
one of the v microcells with no cohabitation. The number 
of ways of assigning a microcell to the first element of a 
macrocell is f , to the second v — 1 etc. Assuming that 
the particles are indistinguishable, the number of ways of 
assigning microcells to all n, particles in a macrocell is 
thus 



1 



ml {v-rii)\ 



(23) 



To obtain the number of microstates corresponding to the 
macrostate {rii} defined by the number of particles rii in 
each macrocell (irrespective of their precise position in the 
cell), we need to take the product of terms such as H23[) 
over all macrocells. Thus, the number of microstates cor- 
responding to the macrostate {rii}, which is proportional 
to the a priori probability of the state {rii}, is 



ni\(v - rij)! ' 



(24) 



This is the Fermi-Dirac statistics which is applied here in 
physical space. As is customary, we define the entropy of 
the state {rii} by 



S{{n i }) = lnW({n i }), 



(25) 



It is convenient here to return to a representation in terms 
of the density in the i-th macrocell 



rii m riiOQ 
vh d v 



(26) 



where we have defined o~q — m/h d , which represents the 
maximum value of p due to the exclusion constraint. Now, 
using the Stirling formula, we have 

\nW({ ni }) ~ V vftnv - 1) - v{— [lnf —\ - 1 



i_£L 

CO 



ln<M 1-^ 

CO 



(27) 



Passing to the continuum limit v — > 0, we obtain the ex- 
pression H21(l of the Fermi-Dirac entropy in physical space. 
In the dilute limit p <C o~o, it reduces to the Boltzmann 
entropy Q. 

The entropy is the correct thermodynamical potential 
for an isolated system for which the energy is conserved 
(microcanonical ensemble). This is not the case for our 
system which is dissipative. The proper description is the 
canonical ensemble and the correct thermodynamical po- 
tential is the free energy F = E — T e fjS constructed with 
the Fermi-Dirac entropy (1221) and the energy (|21|l . The 
equilibrium state in the canonical ensemble is obtained by 
minimizing the free energy at fixed mass. Writing the first 
variations as SF — XSM = 0, we obtain 



SF 
J~p 



A = 0, 



(28) 



which leads to the Fermi-Dirac distribution (|20|l with a, = 
—A/3. Now that the proper thermodynamical potential has 
been derived by a combinatorial analysis, we can introduce 
phcnomcnologically a dynamical model by writing the evo- 
lution of the density as a continuity equation dtp = V • J 
where the current is the gradient of the functional deriva- 
tive of the free energy, i.e. 



dp _ / SF 



(29) 



This formulation ensures that the free energy decreases 
monotonically provided that p > 0. Indeed, 

• f SFdp , f SF^ , 

= - J J -V^dr = - J p(^J^J dr<0. (30) 

Furthermore, a steady state satisfies F = 0, i.e V(8F/Sp) = 
leading to Eq. (J2SJ. Now, using Eqs. (T2TJ and (T22J), we 

have 



Sp 



-Vc- 



D 



X P(l-p/o-o) 



(31) 



To avoid the singularity when p = or p = o"o > w e require 
that n is proportional to p{\ — p/<Jo)- Writing Eq. (|2"9"|l in 
the form 



dt 



SF 

XP(1 ~ p/cr )V — 

Sp 



(32) 
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and using Eq. (|31j) . we obtain Eq. 1(17|1 . This approach to 
construct relaxation equations is equivalent to Onsager's 
linear thermodynamics. Indeed, noting that the potential 



A(r,t) = 



dp 



T eff In 



1 - p/(T 



(33) 



is uniform at equilibrium according to Eq. I|2U|) or (|28[) . 
the linear thermodynamics of Onsager suggests to writing 
the current as 



/iVA(r,t), 



(34) 



which is equivalent to Eq. I|29|> . The same results can 
be obtained by a variational formulation which is related 
to the Maximum Entropy Production Principle (MEPP) 
|15| . The rate of dissipation of free energy is given by 



F = 



S£dp 
Tp~~di 



dr = 



SF. 

Tp 



V • J dr = - 



SF 

J • V— dr. 

dp 



(35) 



We shall determine the optimal current J* which maxi- 
mizes the rate of dissipation of free energy F under the 
constraint J 2 < C(r,t) putting a (physical) bound on |J|. 
The corresponding variational problem can be written 



5F + S I —dr = 0, 

2/j. 



(36) 



where p, is a Lagrange multiplier. Performing the varia- 
tions on J, we obtain J* = pV(SF/Sp) which returns Eq. 



2.3 Kinetic derivation of the model 

As discussed previously, Eq. 114(1 can be viewed as a non- 
linear Fokker-Planck equation where the diffusion coef- 
ficient and the mobility explicitly depend on the local 
concentration of particles. Such generalized Fokker-Planck 
equations can be derived from a kinetic theory, starting 
from the master equation, and assuming that the proba- 
bilities of transition explicitly depend on the occupation 
number (concentration) of the initial and arrival states. 
Below, we briefly summarize and adapt to the present sit- 
uation the approach carried out by Kaniadakis 26 in a 
more general context. 

We introduce a stochastic dynamics by defining the 
probability of transition of a bacteria from position r to 
position r'. Following Kaniadakis we assume the fol- 
lowing form 

7r(r r') = w(r, r - r>[p(r, t)}b[p(r' ,t)}. (37) 

Usual stochastic processes correspond to a(p) — p and 
b(p) = 1: the probability of transition is proportional to 
the density of the initial state and independent on the 
density of the final state. They lead to the Fokker-Planck 
equation as will be shown below. Here, we assume a 
more general dependence on the occupancy in the initial 



and arrival states. This can account for microscopic con- 
straints like close-packing effects that can inhibitate the 
transition. Quite generally, the evolution of the density 
satisfies the master equation 



dp 
dt 



= / [7r(r' — > r) — n(r 



\dr'. 



(38) 



Assuming that the evolution is sufficiently slow, and local, 
such that the dynamics only permits values of r' close to r, 
one can develop the term in brackets in Eq. I|38|l in powers 
of r — r'. Proceeding along the lines of we obtain a 
Fokker-Planck-like equation 



dp 
dt 

with 



and 



d 

dx, 



, \ , , d In k(p) » dp 



j(p) = a(p)b(p), n(p) = 



a (p) 

b{ P y 



Ci(r) = - J yitu(r,y)dy, 



(39) 
(40) 

(41) 
(42) 



The moments Q and Qj are fixed by the Langevin equa- 
tions JSj. Assuming isotropy 



Ci J%i Cij — F)5ij , 

the kinetic equation becomes 



dp 

at 



= V • 



(J + V^) 7 (p)+ 7 (p)^^^Vp 



(43) 



(44) 



Now, according to the Langevin equations JSJl, D is inde- 
pendent on r and J = — x^ c - Thus, we get 



dt 



(91nK(p) 
D 1\P) q V P - X7(p)Vc 



If we define 

Hp) = lip) 



dp 

a in Hp) 
"dp— 



g(p) =Hp), 



the foregoing equation can be written 
dp 



dt 



\7-[Dh(p)Vp- X g(p)Vc], 



(45) 



(46) 



(47) 



and it coincides with the phenomenological equation l|13l) . 
It seems natural to assume that the transition probability 
is proportional to the density of the initial state so that 
a(p) — p. In that case, we obtain an equation of the form 



dp 

at 



V-(D[b(p)-pb'(p)]Vp-xpb(p)Vc) 



(48) 



Note that the coefficients of diffusion and mobility are 
not independent since they are both expressed in terms 
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of b(p). Choosing b(p) — 1, i.e. a probability of transition 
which does not depend on the population of the arrival 
state, leads to the standard Fokker-Planck equation ©. If, 
now, we assume that the transition probability is blocked 
(inhibited) if the concentration of the arrival state is equal 
to (To, then it seems natural to take b(p) = 1 — p/<r - In 
that case, we obtain 

^ = V • {DVp - X p(l ~ p/a )Vc) , (49) 

which coincides with the phenomenological equation (|17|) . 

We can consider a related kinetic model with simi- 
lar thcrmodynamical properties and the same equilibrium 
states. For this model, the dynamical equation reads 



dp 
di 



T, 



If 



1 - p/co 



Vp - pVc 



(50) 



This can be put in the form of a generalized Smoluchowski 
equation |15| : 



dp 
dt 



= V • 



(Vp - pVc) 



(51) 



associated with a barotropic equation of state p(p) = 
— <7oT e ff ln(l — p/ao), where p is an effective "pressure". 
For p <C (To i we recover the "isothermal" equation of state 
p = pT e f f leading to the ordinary Keller-Segel model (J3J) . 

Equations l|49|l and Ij50|l have very similar properties 
and they can be viewed as natural extensions of the Keller- 
Segel model. In Eq. I)49|l the regularization is put in the 
drift term while in Eq. (| 5 f )|> it is put in the diffusion term. 
These two possibilities are considered in JS]. Note finally 
that Eq. JSOl can be obtained from the master equation 
Ij38|l when the transition probabilities are of the form l|37|) 
with a(p) = pj yjl - p/a and b(p) = \/l - p/cr . 



3 Domain wall theory 
3.1 The stationary state 



Introducing the new variables 
c a 



(Aao/2)' 



C 



4fc 2 ' 



we get 



Aip - k 2 ip 



-1 + tanh 



C 



(g-*v) 



(55) 



(56) 



We shall solve this equation perturbatively in a limit of 
high degradation so that k ^> 1 (note that the following 
treatment is exact in ID without approximation). In that 
limit, we can neglect the gradient of concentration except 
in a thin layer of length ~ fc" 1 (domain wall) where the 
concentration changes rapidly. Outside the wall, we obtain 
the algebraic equation 



- 1 + tanh 



C 



\c 



k 2 Tp 



(57) 



In some range of parameters (see below) , this equation de- 
termines two solutions ip± which correspond to two phases 
with uniform concentration. These two phases are con- 
nected by a "wall". For k 3> 1, the interface is very thin 
so that we can neglect the curvature of the wall in a first 
approximation. The wall profile is then determined by the 
one-dimensional differential equation 



" ^ = - 1 + tanh [ c - k ^ 

where £ is a coordinate normal to the interface. 

3.2 The wall equation 

If we set 

(j) = k 2 ip- p/C, T = k£, x= — -1, 



(58) 



we obtain the wall equation 



The stationary solution of the regularized Keller-Segel model 
l|18|l is the Fermi-Dirac-like distribution 



= - tanh(C0) + <t> + x- 



(59) 



(60) 



P 



CO 



1 + e^l 3c+a ' 



(52) 



Eq. (|60f) is similar to the equation of motion for a particle 
in a potential 



This relates, at equilibrium, the bacterial density p to the 
concentration of the chemical c. The chemical is itself pro- 
duced by the bacteria according to Eq. (fTTTfl . Thus, com- 
bining Eqs. H52|l and (|19J) . we obtain a differential equation 
for the concentration c. Using the identity 



U{<j>) = Tin [cosh(0/T)] - ^ + X <t> + U , (61) 

where r plays the role of time and (f> plays the role of 
position. Indeed, it can be rewritten 



1 



1 



1 — tanh 



this mean-field equation can be written 

' a — fic 



Ac - k z c = 



Aero 
2 



1 — tanh 



(53) 



(54) 



dT 2 



-u'(4>). 



(62) 



We have introduced the notation C — 1/T where, as we 
shall see, T plays the role of a temperature (it is further- 
more proportional to the effective temperature T e ff). On 
the other hand, Uq is a constant of integration that will 
be specified later. 



P.H. Chavanis: Phase separation of bacterial colonies in a limit of high degradation. 



7 



Far from the wall, the density is uniform with values 
4>± satisfying 

U\4>±) = 0. (63) 
On the other hand, a first integral of Eq. 1|62|) is 

The condition of solvability is therefore 

U{cj>.) = U(<f>+). (65) 

The only possibility to satisfy the two conditions (|63|l - 
(|65l) simultaneously is that x = 0> i- e M — ft ( m that case, 
E/(0) is symmetric and the above conditions are satisfied 
trivially). Then, the wall profile is completely determined 
by the equations 

^|=-tanh(0/T) + 0, (66) 



T>1 

- 


;t=tc/ / 
/ / / T<Tc 


Iff 









-4 -2 2 4 



Fig. 1. Graphical construction determining the solutions ±u of 
the algebraic equation l|68|l -b as a function of the temperature 
T. 

2 | , , , , , , , , , 1 



U{4>) = T In [coshfo/T)] - y + C/o, (67) 



= ±u, 



u = tanh(w/T). 



(68) 



We note that the algebraic equation has solutions m ^ 
only if (see Fig. ^| 



T < T r = 1. 



(69) 



This is similar to a second order phase transition (see Fig. 
|2J . In that case, the algebraic equation has three solutions 
and ±u with u < 1 (but the solution = tt = is 
unstable). The concentrations of bacteria and chemical in 
the uniform domains are related to the order parameter u 

by 



ft _1_ N 



P± 



(TO 



(lzbtl). 



(70) 



More generally, using Eqs. lj521) - (|S"3l the concentration pro- 
files in the whole space can be expressed in terms of the 
field 4> by 



2k 2 



(1 + 0), p = ^ [1 - tanh^/r)] • (71) 



1.5 
1 

0.5 


-0.5 
-1 
-1.5 



\ T=Tc 



u=0 



T<Tc 



T>Tc 



0.4 



0.8 1.2 

T/Tc 



1.6 



Fig. 2. Evolution of the order parameter u as a function of the 
temperature T. For T > T c the system is homogeneous with 
<j> = u = 0. For T < T c , the uniform solution becomes unstable 
(see Appendix [XJ and two phases <j>± = ±it separated by a 
"domain wall" appear. 



With this convention, the constant appearing in Eq. (|64|l 
is E = 0. Then, we obtain the equation 



2 I* =" £ "«- 



(73) 



3.3 The wall profile 

We determine Uq such that E/(u) = 0. The potential is 
then explicitly given by (see Fig. : 



U{4>) = Tin 



cosh — 
T 



- — - Tin 
2 



cosh 



(f) 



Y' 

(72) 



which determines the wall profile by a simple integration 
(see Fig. 



fix 



o v/-2(7(x) 



(74) 



For r - 
we set 



+oo, 
= u - 



) — > u. To get the asymptotic behaviour, 
with 6 <^ 1 and we linearize the wall 
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-0.5 - 



T=0.7 


/ 
// 

// 

// 

/ 

ii 

1 


/ 
// 
// 
ff 

// 

-u 





-10 



10 



Fig. 3. The potential U(<j>) of the equivalent mechanical prob- 
lem for different values of the temperature. 



Fig. 4. Plot of the field (f> accross the wall for T = 0.7. The solid 
line corresponds to the exact solution of Eq. 173H obtained nu- 
merically and the dashed line corresponds to the approximate 
expression HOUt . 



equation (|73")l . This yields 



cW 



(75) 



with U"(u) = i(l - u 2 ) - 1 < for T < T c . The wall 
connects the uniform phase exponentially rapidly Thus, 
we can write 



<f> = u- A{u)e- 2T/L(u) , (t^+oo) 



(76) 



where the typical width of the wall (expressed in units of 
fc- 1 ) is 



0.4 



0.3 



0.2 



0.1 - 



T=0.7 



L(u] 



We introduce the concentration gradient (see Fig. EJ) 



(77) 




v(r) = f r = V^UW). 



(78) 



Using Eq. J72J, we find that the maximum value of the 
concentration gradient, corresponding to (f> = 0, is given 

by 



Umax (^0 



2T In 



cosh 



(79) 



Fig. 5. Plot of the concentration gradient v(t) accross the wall 
for T = 0.7. The solid line corresponds to the exact solution of 
Eq. I73H obtained numerically and the dashed line corresponds 
to the approximate expression 11031 - 11051 . 



3.4 The limit T -> T c 

For T — > T c , u — > and (f> <C 1. In that case, we can 
expand the potential to order (j> to obtain 



Finally, the energy of the wall, or surface tension, is 



U^)= l -{\-T)4> 2 -^ + U,. 



a(u) 



7h 



dT 



+ u 



y/-2U{<j))d<j). (80) 



The functions L, a and v max are plotted in Figs. |H| El 
and|Slas a function of the temperature T, together with 
their asymptotic expressions computed in the following 
sections. The concentration profiles of bacteria and of the 
secreted chemical, given by Eq. l(7T|l . are plotted in Fig. El 



The wall equation becomes 
d 2 6 



dT 2 



(l-7>+-0 3 . 



(81) 



(82) 



The uniform solutions are tfi — and <fi 2 — 3(1— T) yielding 



(83) 
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\ X X 




















\ 
\ 



0.2 



0.4 0.6 
T/Tc 



0.25 



0.5 
T/Tc 



0.75 



Fig. 6. Evolution of the typical width of the wall as a func- 
tion of the temperature. The dashed lines correspond to the 
asymptotic expressions for T — > and T — > T c . 

1.5 i ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 



Fig. 8. Evolution of the energy of the wall (surface tension) 
as a function of the temperature. The dashed lines correspond 
to the asymptotic expressions for T — > and T — > T c . 




0.5 - 




Fig. 7. Evolution of the maximum concentration gradient as 
a function of the temperature. The dashed lines correspond to 
the asymptotic expressions for T — > and T — > T c . 



We can now re-express the potential in the form 

The wall profile is given by 

d<j> t 



Fig. 9. Concentration profiles of bacteria p (in units of cro) and 
of the secreted chemical c (in units of Acro/2fc 2 ) for T = 0.7. 



The width of the wall diverges at the critical point with 
the exponent —1/2. The wall profile can be rewritten 



yielding explicitly 



u 2 - ^ ^6' 



. . UT 

= u tanh — = 
.\/6 



(84) 



(85) 



(86) 



<j> — u(T) tanh 



|L(T) 

and the concentration gradient is 



I \ U (T) i -2 
V[T) = ; cosh 



L(T) 



L(T) 



(89) 



The typical width of the wall, as defined by Eq. I)77|) . is 
\/6 nz.„ ^-1/2 



The maximum value of the concentration gradient is given 

by 



L = — = V2(l -T)- 1 ' 
u 



(87) 



1/2 



(i-n 



(90) 
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and it tends to zero with the exponent +1 at the critical 
temperature. Finally, the surface tension is given by 



3/2 



[2(1 - T)] 3/2 , 



(91) 



and it vanishes at the critical point with the exponent 
+3/2. These are the same scalings as in the Cahn-Hilliard 
theory [27|. 



3.5 The limit T -> 

Setting x — u/T, Eq. (|68|) -b can be rewritten Tx = tanh(x). 
For T — > 0, x ~ 1/T — > +oo and u — * 1. More precisely, 
considering the behavior of tanh(x) for a; — > +oo, we get 

u ~ 1 - 2e _2/T . (92) 
The potential can be rewritten (for > 0) 



(93) 



The typical width of the wall is 



L = 2 (l + | e - 2 / r ) = 2 - (1 - «) In (i^) , (94) 
and the maximum value of the concentration gradient is 

(95) 



, , 2 In 2 
1 - ln2 T = 1 j-, r . 



Finally, the surface tension is given by (see Appendix [BJ) 



a = 1 T = 1 - 



12- 61n(^)' 
For T = 0, we have u = 1 and 

4> 2 1 l . , 
u(4) = <l>- =- 5 W-i) 2 . 

The wall profile is solution of 

( /c6 



fir 



1 



leading to 



= l-e" r (t>0). 



(96) 



(97) 



(98) 



(99) 



3.6 Simple approximation for T > 1/2 

If we are relatively close to the critical temperature, we 
can propose a simple approximation of the wall profile in 
the form 



where u and L are given by the exact expressions 

u = tanh (|) , (101) 

2 



L 



i-A(i 



(102) 



This Ansatz becomes exact when T — ► T c and it provides 
a fair approximation for smaller temperatures (typically 
T > 1/2), see Fig. |U 

The concentration gradient obtained from Eq. (|10Ufl is 
given by 



<p = u(T) tanh 



L(T) 



(100) 



v{t) 



-XT) 



cosh 2 '{t/L{T)Y 
with the maximum value 



u u 1 . 

- = -\/l (1 - m 2 )- 

L 2 V T 



(103) 



(104) 



However, it is more relevant to take for v max the exact 
value (f?9")) . i.e. 



2Tln 



cosh ( 

T 



(105) 



Typically, Eq. i|103|) with l|104|l gives a better agreement 
with the exact solution in the tail of profile while Eq. i|103[l 
with l|104(l gives a better agreement in the core of the 
profile, see Fig. [3J Finally, the surface tension calculated 
with Eq. (|103|l is given by a 



4 2 t 
3 v max^ i ' 



3.7 Simple approximation for T < 1/2 

For sufficiently small temperatures, we can propose a sim- 
ple approximation of the wall profile in the form 



4> = u(T) 



1 



,-2t/L{u) 



(r>0), 



(106) 



where u and L are given by the exact expressions <|101ll 
and 1|102|) . The concentration gradient obtained from Eq. 

CEDE) is 



u(t) 



XT)e 



-2\t\/L(T) 



(107) 



with the maximum value v max = 2-^. As before, it may 
be more relevant to use the exact value (|105|) . The surface 
tension calculated with Eq. (|107fl is given by a = \v^ lax L 
but this approximation yields, unfortunatelly, an asymp- 
totic behaviour for T — > different from the exact result 

The Ansatz <|100[) - (110611 have a self-similar structure 
as a function of the temperature. Indeed, the functions 
4>(t)/u(T) and v(T)/v max (T) vs t/L(T) have an invari- 
ant profile. The exact solution of the differential equation 
(|66|l has not an exact self-similar structure as shown in 
Figs. 1101 and 1111 but the region between the envelopes is 
relatively thin so that Eqs. H100|) - (|106fl can be useful ap- 
proximations for T — > 1 and T — > respectively. The 
profiles of concentration and concentration gradient (in 
non-scaled variables) are plotted in Figs. 1121 and 1131 
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1-exp(-2x) 


a ^ tanh(x) 


Jm 






-0.5 = 



-5 5 T 

x/L 

Fig. 12. Concentration as a function of the distance r for 
Fig. 10. Scaled concentration <f>/u(T) as a function of the different values of the temperature T = 0.1, 0.2, 0.3, 0.4, 0.5, 
scaled distance r/L(T) for different values of the temperature rj.6 q 7 q g anc [ q g 
between T = and T = T c = 1. In terms of the scaled vari- 
ables, the exact concentration profile is bounded by the solu- 1 
tions 1 - cxp(-2:r) for T = and tanh(x) for T = T c . 



> 0.5 - 









^ 1/cosh 2 (x) 


If 1 

#1 






i - 




-5 





T/L 



Fig. 11. Scaled concentration gradient v/v ma x(T) as a func- 
tion of the scaled distance r/L(T) for different values of the 
temperature between T = and T = T c = 1. In terms of the 
scaled variables, the exact profile of concentration gradient is 
bounded by the solutions exp(— 2x) for T = and 1/ cosh 2 (a;) 
for T — T c . 



3.8 Match asymptotics 

For small values of the temperature, one can propose an- 
other approximation of the profiles of concentration and 
concentration gradient by using match asymptotics. 



3.8.1 Concentration profile <p(r) 

The exact asymptotic behaviors of the concentration pro- 
file are given by 



Fig. 13. Concentration gradient v as a function of the distance 
r for different values of the temperature T = 0.1, 0.2, 0.3, 0.4, 
0.5, 0.6, 0.7, 0.8 and 0.9. 



(r) = u - Ae- 2r l L (t -» +00), 



(109) 



where u, L and v max are known functions of the temper- 
ature. We match these two behaviors at a point x where 
their values and the values of their first derivative coin- 
cide. This yields 



Til — II — Ap~ 2 ' X / L 

v max — c - 

From these relations, we obtain 

u L 

Umax 2 



X 



</>(t) ~ v max T+ ... (r->0), 



(108) 



(110) 

(111) 

(112) 
(113) 
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0.8 



0.6 - 



0.4 



0.2 



T=0.3 




Fig. 14. Concentration profile for T = 0.3: exact (solid), Fig. 15. Concentration gradient for T = 0.3: exact (solid), 
approximate self-similar (dotted), match asymptotics (long- approximate self-similar (dotted), match asymptotics (long- 



dashed) . 



dashed) . 



For T 



0, we explicitly find that 
x = (In 2)T, 



A=l- (ln2) 2 T 2 



(114) 



(115) 



The exact concentration profile is plotted in Fig. El f° r 
T = 0.3, and compared with the approximate self-similar 
expression i JTUffl and the expression (|108[l - i|109|) obtained 
by match asymptotics. 



3.8.2 Concentration gradient v(t) 

The exact asymptotic behaviors of the concentration gra- 
dient are given by 



v(t) ~ v n 



(r-0), 



v(t) = 2ye- 2T ^ L (t -> +oo). 
Li 



(116) 



(117) 



We match these two behaviors at a point x where their 
values and the values of their first derivative coincide. Af- 
ter simplification, we obtain 



—L 



L- 



8T 
1-T 



1 



For T — > 0, we explicitly find that 
£ = T, 



(118) 
(119) 

(120) 



A = l + (l-ln2)T. 



(121) 



The exact concentration gradient is plotted in Fig. I15lfor 
T = 0.3, and compared with the approximate self-similar 
expression (|107[l - (|105[) and the expression (|ll(j[) - l|117|) ob- 
tained by match asymptotics. 



3.9 The curvature radius 

To next order in the expansion in fc _1 -C 1, we must ac- 
count for the curvature of the interface. Close to the in- 
terface, VV> = f^n where n is a unit vector normal to the 

wall. Introducing the curvature = Vn (r is the curva- 
ture radius), we get Atp 
Eq. <|56(l becomes at first order 



d 2 ip/d^ 2 +r- 1 dtp/d^. Therefore, 



dr 1 



1 d(j) 
kr dr 



-u'(4)+xi/k, 



(122) 



where we have written X — Xo + k~ 1 xi + ■■■ and used the 
fact that xo = at leading order (see Sec. 13.3(1 . Multiply- 
ing Eq. I|122|) by dtfi/dr and integrating along the wall we 
obtain the relation 



2xm 

a(u) ' 



(123) 



which shows that the radius of curvature r is constant. 
Therefore, the shape of the interface is either a circle (lead- 
ing to a spot) or a straight line (leading to a stripe). Equa- 
tion (|123fl is similar to Laplace's law relating the curvature 
radius of bubbles to the surface tension and to the differ- 
ence of pressure between the interface. 



3.10 The free energy 

The previous results can be recovered by minimizing the 
'free energy' l)21 [l -l)22 [l . This has been discussed by Bouchet 
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& Sommeria |24j in the context of jovian vortices where 
the problem is similar (see Sec.0J. Therefore, their study 
can be directly applied to the present situation and we 
shall only give the main steps of the analysis. 

To leading order in fc _1 — > 0, the system consists of 
two phases with uniform density p±\ and size A±\ (we 
call A = A + + A- the total size of the domain). Setting 
p = ^(1 + (/)) and c = + <p), the free energy F = 

E — TeffS + aT e ffM where E and S are given by Eqs. 
J2U-(|22l and M is the total mass can be written 



F = A 1 f(p 1 ) + (A-A 1 )f(p_ 1 ), (124) 



with 



f= T^ao { _ C(f)2 + (a _ 2C)(f) 

+ (1 + (j)) ln(l + <f)) + (1 - 0) ln(l - <)>)}. (125) 

The optimal values of p±\ and A±\ are obtained by mini- 
mizing the free energy 1124(1 . The variations on p±\ imply 
that 



f(p±i) = 0, f"(p±i) > 0. 
and the variations on A\ imply that 

/(Pi) = /(p-i)- 



(126) 



(127) 



This relation expresses the equality of the free energy of 
the two phases. The only possibility to satisfy the relations 
<|126|) - (I127[I simultaneously is to have a = 2C so that /(</>) 
is an odd function. This is equivalent to the solvability 
condition 1(65(1 leading to X = otjlC —1 = 0. Then, it 
is straightforward to check that f'(p) — implies that 
tj> = ±it where u is given by 



Cu = - In 
2 



1 



= tanh (u). 



(128) 



This returns the relation l(68(l -b. 

To first order in fc _1 , we need to determine the contri- 
bution of the free energy contained in the wall (interface) . 
The free energy per unit length is given by 



F 



w 



[h(p(T))-h(p ±1 )]dr, 



(129) 



where h(p) is the density of free energy. Using Eqs. ((21|) . 
1(221) . 1(71(1 and ((73(1 . we obtain after simplification 



F 



IV 



4fc 3 



h{4>) - ~h(<p 



±i, 



G?T, 



with 



/i((/>) = 0(tanh(C0) - (j>). 



(130) 



(131) 



Using h((j>±i) — h(u) — and tanh(C</>) > (see Fig. 

we find that Fyy > 0. Therefore, minimizing the free 
energy amounts to minimizing the length of the interface 
at fixed area. This gives either a circle or a straight line and 



this returns the fact that the curvature radius is constant 
(see Sec. 13.9(1 . This argument also shows that it is more 
profitable to form, at equilibrium, a single "bubble" of size 
A rather than several "droplets" of smaller size. However, 
a configuration with several "droplets" can exist as a non- 
equilibrium solution of the regularized Keller-Segel model 
((18(1 - 1(19(1 . These droplets will evolve in time and merge 
together to finally form a single "bubble" (spot or stripe) . 
This is similar to a coarsening process in spin systems or 
to the aggregation of vortices in 2D decaying turbulence 



3.11 Bifurcations: spots and stripes 

We shall work in a square domain with total size A. We 
consider periodic boundary conditions in order to avoid 
boundary effects. The equilibrium state consists in two 
phases with uniform density (p+,A + ) and (p_,v4_) with 
A + + A- = A. We introduce the parameter 



B = 1 



2M 
vqA' 



(132) 



Since < M < ctqA, the parameter B takes values be- 
tween — 1 and +1. Writing M = A + p + + A_p_ and using 
Eq. l|70"|) . the area of the two phases at equilibrium are 
given by 



.4 
2" 



IT 



B 



(133) 



They are determined by the parameter B (fixed by the 
total mass) and by the parameter u (fixed by the temper- 
ature). Since < A± < A, we have the inequalities 

\B\<u<l, -1<B<1. (134) 

Since the curvature radius is constant, we have three pos- 
sible configurations: (i) a circular domain (spot) of phase 
+ surrou nded by phase — : the length of the interface is 
2WttA+. (ii) a circular domain (spot) of phase — sur- 
rounded by phase +: in that case, the length of the in- 
terface is 2y>'ir(A — A + ). (iii) a stripe of phase ± and a 
stripe of phase =F: the length of the interface is 2^/A. The 
configuration selected at equilibrium is the one with the 
smallest interfacial length (see Fig. 116(1 . 

(i) The spot of phase + surrounded by phase — will 
be selected if A + < A/ir, i.e. u < irB/(n - 2). This is 
possible only for B > 0. In terms of the temperature, 
this corresponds to T > To where Tq(B) is determined by 
Mo = tanh(uo/To) with Uo = \B\tt/(tt — 2). 

(ii) The circular domain of phase — surrounded by 
phase + will be selected if A + > (1 — l/n)A, i.e. u < 
—irB/ (tt — 2). This is possible only for B < 0. In terms of 
the temperature, this corresponds to T > To- 

(iii) The stripes will be selected for T < Tq. This is 
possible only for \B\ < B c = (ir — 2)/ir (i.e. uo < 1). For 
B > 0, the stripe + has the smallest area (A+ < A_) and 
this is the opposite for B < 0. 
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Fig. 16. Interfacial length of the spot (circle) as a function of 
the temperature for different values of the parameter B (the 
area of the domain has been normalized to A — 1). The spot is 
selected if its length is smaller than 2 and the stripe is selected 
in the other case. If B > B c = (ir — 2)/7r, the spot is always 
selected. If B < B c , the stripe is selected for T < To(B) and 
the spot is selected for T > Tq(B). 
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Fig. 17. Phase diagram of the regularized chemotactic model 
p8(l - l|19| l showing the bifurcation between 'spots' and 'stripes' 
as a function of the control parameters (B,T). In a bounded 
domain, the solutions exist only for T < T ma x(B). We have 
drown the line of transition To(B) between the two structures. 
Finally, the dashed line corresponds to the domain of validity of 
our perturbative expansion; it has been plotted for Ak 2 — 1000. 



Finally, we note that the condition \B\ < u implies 
that T < T max (B) < T c where T max (B) is determined by 

ax J • 

The phase diagram is represented in Fig. El This is the 
counterpart to the diagram obtained by Bouchet & Som- 
meria for jovian vortices. The 'spots' are the equiva- 
lent of the 'vortices' and the 'stripes' are the equivalent of 
the 'jets'. The main difference (beyond the context and the 
interpretation of the solutions) is that the control parame- 
ter in our case is the effective temperature T ( "canonical" 
situation) while their control parameter is the energy E 
(microcanonical situation). 

Finally, concerning the domain of validity, our approach 
assumes that the domain size is larger than the interfacial 
width L(u)/k given by O so that A± > (2L/k) 2 . This 
is satisfied on the left of the dashed line in Fig. El corre- 
sponding to 

l^.(.-<^). 

4 Analogy with Jupiter's great red spot 
4.1 The physical context 

Jupiter's Great Red Spot (GRS) is probably the most fa- 
mous example of vortex structures found in planetary at- 
mospheres. The presence of this spot was first reported 
by Robert Hooke in 1664 in the first issue of the Philo- 
sophical Transaction of the Royal Society. The GRS is 
a large "eye" that dominates the southern hemisphere 
of the planet. It is an oval-shaped anticyclone with size 



26000km by 13000km. Its breadth is about one hundred 
times its height so it can be considered approximately 
two-dimensional. It resides in a zonal shear at latitude 
23S where the velocity changes sign, is elongated along 
the shear zone and is of the same sign as the background 
shear (these are relatively general rules observed for other 
vortices). It stays at the same latitude but slowly drifts 
in longitude. Morphologically, the GRS has a striking an- 
nular structure with a quiet center surrounded by a thin 
intense jet. The concentration of winds in an annulus is 
consistent with the fact that the GRS is much larger than 
the atmosphere's radius of deformation (the jet's width 
scales with the Rossby radius of deformation). 

Some authors have tried to describe this spot as a soli- 
ton solution of a modified KdV equation but this interpre- 
tation as a weakly nonlinear structure faces several draw- 
backs: the jet structure is not reproduced and the soli- 
ton theory predicts an interpenetration without change of 
structure while the main interaction between vortices is 
a merging process. In contrast, Marcus |55] argues that 
a permanent vortex can coexist with turbulence and that 
most of the properties of the Jovian vortices can be eas- 
ily explained and understood with the quasi-geostrophic 
(QG) theory. In addition, he emphasizes the jet structure 
of the GRS and shows that an annular jet is the natural 
structure of a vortex with a uniform potential vorticity in- 
side and outside the spot. A physical justification of this 
construction was given by Sommeria et al. |23| in terms 
of statistical mechanics. A statistical theory of 2D tur- 
bulence has been developed by Miller 19 and Robert & 
Sommeria 20 following and generalizing the pioneering 
work of Onsager |30j on point vortices. This theory is able 
to account for the typical vortices (monopoles, dipoles, 
tripoles,...) observed in two-dimensional flows |H3- I n the 
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QG model, and in the limit of small Rossby radius, PV 
mixing (entropic effects) with constraints on the energy 
leads to an equilibrium state that consists of two phases 
with uniform PV in contact separated by a strong jet. This 
precisely account for the morphology of the GRS. This 
model has been further developed by Bouchet & Somme- 
ria |24| with quantitiative applications to Jovian vortices. 
Another model has been proposed by Turkington et al. 
|32| . It predicts the emergence of a vortex solution at the 
correct latitude but it does not reproduce the annular jet 
structure of the GRS 33 . 



4.2 Statistical mechanics of the quasi-geostrophic 
equations 

The quasi-geostrophic equations appropriate to the dy- 
namics of geophysical flows [22] can be written: 



q 



dq 
dt 

± 
R 2 



u- V<? = 0, 



(136) 



Rh(y), u=-zxVf (137) 



Here, q is the potential vorticity (PV) and ip the stream 
function (z is a unit vector normal to the two-dimensional 
flow). We have assumed that the topography Rh(y) scales 
with the Rossby radius R. The Q.G. equations admit an 
infinite number of stationary solutions specified by any 
relationship q = f(ip). For given initial conditions, the 
statistical theory selects the most probable state consis- 
tent with the constraints imposed by the dynamics. It is 
obtained by maximizing a mixing entropy at fixed energy 
and Casimir constraints 19 20,21. 

Let us consider the situation where the fine-grained 
PV q takes only two values {a_i,ai}. It is convenient to 
set (ai — a_i)/2 = 1 and (ai+a_i)/2 = B. We also choose 
the Gauge condition on ip such that (q) = ~~~J. Then, 
the total area occupied by level a\ is A = (1 — B)/2 (the 
total area of the domain is unity). In this two-levels case 
|19U20U21| . the mixing entropy is given by 



S = - J \pbxp + (1 - p) ln(l - p)]dr, 



(138) 



where p(r) is the local distribution of level a\ . The coarse- 
grained PV is q — pai + (1 — p)a~i. The extremization of 
Ijl38(l at fixed energy 



dr, (139) 



and total area A = J p(r)dr leads to a q — t/j relation of 
the form 



Ctp 

q = B- tanh ( a - — ) , 



(140) 



where a and C are Lagrange multipliers introduced in the 
variational problem SS — 2aSA + jpSE = 0. 



Robert & Sommeria |34| have proposed a parameteri- 
zation of 2D turbulence in the form of a relaxation equa- 
tion for the coarse-grained P V q(r, t) . This parameteriza- 
tion is based on a Maximum Entropy Production Principle 
(MEPP). The diffusion current (due to turbulent mixing) 
is assumed to maximize the rate of entropy production 
S while conserving all the constraints imposed by the dy- 
namics. In the two-levels case, this yields a system of equa- 
tions of the form 

dq _ 

V ■ [D (Vg + (3(t)(a x -?)(?- a_i)VV0] . (141) 



J DVq ■ S7ip dr 



J D(ai -g)(g-a-i)(V*/0 2 dr' 



q = -A^ + ^-Rh(y). 



(142) 



(143) 



Interestingly, these drift-diffusion equations are similar to 
the regularized chemotactic model Ijl8|l - (|19|l . In this anal- 
ogy, the coarse-grained PV q plays the role of the bacterial 
concentration and the stream function ip the role of the 
chemical c. This analogy was noted in ~~|. An important 
difference, however, is that in 2D turbulence the energy is 
conserved so that the inverse temperature (3(t) evolves in 
time. By contrast, in the chemotactic problem, we are in 
a situation where the 'effective temperature' T e ff = 1/(3 
is fixed. Therefore, 2D turbulence corresponds to a micro- 
canonical situation (where we maximize the entropy <S* at 
fixed energy E) while chemotaxis corresponds to a canon- 
ical situation (where we minimize an effective free energy 
F = E — T e ffS). These two situations are described in 
|15j. For long-range interactions (k = or R — > +oo), the 
'ensembles' are generically inequivalent. However, in the 
limit that we consider here, corresponding to a short-range 
interaction, they become equivalent, i.e. the (3(E) curve is 
univalued (see Fig. [T5)l . 

The equilibrium problem (|140|) - (|143fl has been studied 
in the limit of small Rossby radius R — * 0. The original 
idea dates back to Sommeria et al. who understood 
that, in this limit, the solution is made of two uniform 
PV regions separated by a strong jet. This model has 
been developed quantitatively by Bouchet & Sommeria 
[2~4] with comparison to jovian data. Our study of the reg- 
ularized chemotactic problem has been directly insiped 
by these works. In complement, we have provided in Sec. 
[3] a more thorough study of the wall equation with use- 
ful asymptotic expansions and analytical approximations. 
Due to the analogy between the two problems, these re- 
sults can also be relevant to describe the jet structure of 
jovian vortices, like GRS. In this respect, we recall that 
the GRS corresponds to a typical parameter u in the range 
0.92 < u < 1 .24 so that the limit u -> 1 or T -> of our 
study fSecs. l3~31l3.7l and l3.8fl is particularly interesting in 
that respect. To strenghten the comparison between the 
two problems (chemotaxis and jovian vortices), we briefly 
recall the main lines of the study of Bouchet & Sommeria 
|24j and provide some complementary discussion. 
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4.3 Domain wall theory of Jupiter's great red spot 

4.3.1 The jet equation 

Combining Eqs. (|143[) and (|14Ufl . we find that the stream- 
function satisfies the meanfield equation 

-Ai/j + ^- Rh{y) = B- tanh (a - . (144) 

We shall solve this equation perturbatively as an expan- 
sion in powers of R. We only give the main lines and refer 
to Bouchet & Sommeria [22] f° r more details and devel- 
opements. To first order, we obtain the jet equation 

d 2 v id* i/> ( Ci> 

(145) 

where r is the curvature radius. Introducing the notations 

ip a 



R 2 C" 



(146) 



we get 

-t4 + -? + Rh{v) = - tanh(C<£) + cj) + % - B. 

(147) 

To leading order in R <C 1, the foregoing equation reduces 
to 



d 2 ^ 
dr 2 



tanh(C0) + (f> - 



C 



B. 



(148) 



The condition of solvability implies that ao = CB so that 
the jet equation takes the form 



d 2 ^ 
dr 2 



•tanh(C^) + <f> = -U'((j>). 



(149) 



This is the same equation as Eq. I|66[l with C = 1/T. 
In the present context, <fi is related to the streamfunction 
and v — d(j>/dr to the jet velocity. Therefore, the figures 
representing the profile of concentration gradient in Sec. 
El give the jet velocity profile in the present context. The 
PV and streamfunction in the two phases are 



q± = 

Their area is given by 
A± = 



R 2 



B±u. 



2 I u 



(150) 



(151) 



which is similar to Eq. I|133|) . Finally, in the present con- 
text, the parameter u is determined by the energy accord- 
ing to the relation 



E = X -R 2 (u 2 -B 2 ). 



(152) 







B=B 

c 




• 




• — * 


B=0.3 




B=0.1 B=0.2 







0.2 



0.4 0.6 
2E/R 2 +B 2 



0.8 



Fig. 18. Caloric curve giving the inverse temperature C as 
a function of the energy E for jovian vortices in the limit of 
small radius of deformation R — » 0. This curve is obtained from 
Eq. I|152|l - H68|l in the absence of topography. In terms of the 
variable 2E/R 2 +B 2 , this curve is independent on B. We have 
indicated by a 'bullet' the point of bifurcation (corresponding 
to u = nB/(n — 2)) between a vortex (spot, left) and a straight 
jet (stripe, left), for different values of B: 0.1, 0.2 and 0.3. 



4.3.2 The underlying shear 

To first order in R, the jet equation can be written 

^ + -^+Rh(y) = -U'(<l>)+Rai, (153) 

where we have written a — ao + Ra\ + ... and used ao = 
CB. Far from the jet, we have 



Rh(y) = -U'(4>) +Rat 



(154) 



Writing = u + R5<f), we get 8(f) = (a>i — h(y))/U"(u) so 
that the velocity v s hear — 8<p'(y)e x is given by 



Vshe 



h'(y) 

-U"{u) l-C(l-u 2 )' 



h'(y) 



(155) 



Now, the analysis of the jet equation in Sec. l3.3l shows that 
the function appearing in the denominator of Eq. 11551) is 
related to the jet width (|77ll so that we can write 



Vshear{y) = (y) L ( u ) 2 ■ 



4.3.3 The curvature-topography relation 



(156) 



Finally, multiplying Eq. (|153|l by d<p/dr and integrating 
accross the jet, we obtain 



e(u) 



u(h(y) - qi), 



(157) 
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which relates the radius of curvature r to the underlying 
topography h(y). For a given topography, Eq. Ijl57(l deter- 
mines the form of the jet. This problem has been studied 
in detail in Bouchet & Sommeria in the case of a 
quadratic topography. As another example, we consider 
here the inverse problem: given the form of the jet, find 
the corresponding topography. We consider the case of an 
elliptical vortex (see Fig. I19f) because this is a relatively 
good representation of Jupiter's great red spot and we can 
obtain analytical results in that case. 

The topography leading to an elliptic vortex has the 
form 



H 



i + (y/LY 



3/2 : 



(158) 



where L and H are typical horizontal and vertical lenght 
scales. Assuming that this relation holds for \y\ — > +oo, 
we must take ol\ = in Eq. i|157|) to have a vanishing cur- 
vature at infinity where h — > 0. Therefore, the curvature 
radius of the vortex is given by 



i + (y/LY 



-3/2 



(159) 



where we have denned \ — e(u)/uHL. This is the equa- 
tion of an ellipse with major and minor semi-axis: 



2 ' 



(160) 



These relations assume that x < 1- Now, in the case of 
GRS, the aspect ratio 



a 
b 



1 



(161) 



is close to 2 leading to x = v3/2. The major and minor 
semi-axis are then given by a — 2\^3L and b — V3L. 

Consider now the limit \ ~^ 0. In that case, the width 
of the vortex is small with respect to the horizontal topo- 
graphic length (b <C L) and we can make the quadratic 
approximation 



%) 



rr 3fl 2 
H rrU 



(162) 



This is similar to the situation considered by Bouchet & 
Sommeria |24) . Their parameter d is related to our param- 
eter x by d = (3/2)x 2 - For a quadratic topography with 
(x, d) — > 0, our study shows that the vortex shape is an 
ellipse whose major and minor semi- axis are given by Eqs. 
I|16(J|) . In particular, the aspect ratio behaves like 



1 



(d-Q). 



(163) 



This formula is not applicable for a vortex with aspect 
ratio ~ 2 like GRS. Bouchet & Sommeria |2U assume a 
quadratic topography and solve the curvature-topography 



>• 




Fig. 19. Elliptic vortex with aspect ratio of 2 above a topog- 
raphy of the form l|158|l . 



relation numerically. For large values of d, close to its max- 
imal value d m ax = 4/9, the vortex is not an ellipse. Alter- 
natively, if we assume a topography of the form (|158f) we 
find an elliptic vortex for all the values of x < 1. There- 
fore, the form of the vortex is relatively sensitive to the 
underlying topography. 



5 Conclusion 

In this paper, we have studied the equilibrium states of 
a regularized version of the Keller-Segel model describ- 
ing the chemotactic aggregation of bacterial colonies. The 
regularization is justified physically in order to avoid the 
formation of singularities (Dirac peaks) during the dynam- 
ics and obtain smooth density profiles (clumps) instead. 
This regularization accounts for finite size effects and close 
packing effects. In that case, an equilibrium state exists for 
any value of the control parameter (this is similar to con- 
sidering a gas of self-gravitating fermions in astrophysics 
to avoid complete gravitational collapse |SS1)- We have 
studied furthermore a limit of high degradation k — > +oo. 
In previous works f5TCT7H5ll^lllUlllllll2| . the opposite limit 
k = (no degradation) was considered instead. The inter- 
mediate case of a finite degradation must be studied nu- 
merically (in preparation). However, the asymptotic limit 
k — * +oo allows us to obtain analytical results that permit 
to have a clear picture of the bifurcation diagram (between 
spots and stripes) in parameter space. Furthermore, our 
approach is exact in one dimension, for any value of the 
degradation rate k. Our study shows that the physics of 
the problem is sensibly different whether k = or k ^ 0. 

We have also discussed the analogy between the orga- 
nization of bacteria (in stripes and spots) in the chemo- 
tactic problem and the organization of two-dimensional 
turbulent flows (in jets and vortices) in the jovian atmo- 
sphere. These apparently completely different systems are 
however described by relatively similar equations so that 
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an interesting analogy can be developed between the two. 
In this analogy, the jet structure of Jupiter's great red spot 
can be seen as a 'domain wall' that is similar to the inter- 
face separating two phases in contact, as in the biological 
problem. 



A Stability analysis 

We study the linear dynamical stability of an infinite and 
homogeneous solution of the regularized Keller-Segel model 
<|18|) - (|19[) . The unperturbed solution is such that 



k 2 c = Xp. 



(164) 



Linearizing Eqs. lfT^|l - lfT^|l around this steady state and 
writing the perturbation as 5p = <5/5e lqr e°'*, 5c = <5ce lqr e <Tt , 
we obtain 

XP (l - P /<y a )q 2 5c - (Dq 2 + a)6p = 0, (165) 




Fig. 20. Growth rate of the perturbation as a function of the 
wavenumber for T = 1/2 < T c . 



{q 2 + k 2 )5c - XSp = 0. 



( 166 ) B Surface tension for T -> 



This system of equations admits non-trivial solutions only 
if the determinant is zero yielding the dispersion relation 



n=q2{ X \p(l-p/a a ) _ D 



q 2 + k 2 



(167) 



The system is unstable if a > and stable otherwise. A 
necessary condition of instability is 



\p(l - p/a ) - k 2 > 0. 



(168) 



If this condition is fulfilled, the unstable wavenumbers are 

(169) 



q 2 < nAp(l - poo) - k 2 = q 2 max . 



D 

For u — 0, i.e. p = co/2, we find that the instability 
criterion l|168|) corresponds to 



T < T r = 1. 



(170) 



Therefore, the uniform phase u = is stable for T > T c 
and unstable for T <T C where it is replaced by a 'stripe' 
or a 'spot' (see Fig- • The unstable wavenumbers are 



q 2 <(C -l)k 2 = q 2 n 



(171) 



where we recall that C = 1/T. The growth rate (see Fig. 
120(1 can be written 



The maximum growth rate is obtained for 
q 2 = k 2 (VC-l), 

and its value is 

(t * = Dk 2 (VC- l) 2 . 



(172) 



Using Eqs. (|80|l and (|72|l . the surface tension can be writ- 
ten 

a = 2 1 (<P 2 - 2Tln [cosh(0/T)] 

-u 2 + 2T In [cosh(u/T)]) 1/2 d(t). (175) 

Using 

In [cosh(0/T)] = (j)/T - In 2 + ln(l + e" 2 ^), (176) 
and considering the limit T — * 0, we obtain 

a = 2 f J(l- 4>) 2 - 2Tln(l - e - 2 / T + e~ 2 ^l T )d(^. 
Jo 



(177) 



Setting x = 1 — <j>, this can be rewritten 



(7 = 2^ xyjl-^ln[l + e- 2 / T (e 2x / T - 1)] dx. 



For T — * 0, we obtain 



^ 2 /o X { l ~V^ 



l + e- 2 / T {e 2x ' T ~\) 



(178) 

dx. 
(179) 



(173) Setting y = 2x/T, we find that 
(174) 



2/T 

cr = 1 - 2T / In 



l + e - 2 / T (e^-l 



dy 

y ' 



(180) 
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Setting x — 2/T - y, we get 25. W. Jager, S. Luckhaus, Trans. Amer. Math. Soc. 329, 819 

(1992). 



2/T 



1-2T / In 



2/T/ 2/T-x 



For T — > 0, we finally obtain 



dx 26. G. Kaniadakis, Physica A 296, 405 (2001). 

2 / T _ ■ 27. J.W. Cahn and J.E. Hilliard, J. Chem. Phys. 28, 258 

X (1958). 

( 181 ) 28. C. Sire and P.H. Chavanis, Phys. Rev. E 61, 6644 (2000). 

29. P.S. Marcus, Annu. Rev. Astron. Astrophys. 31, 523 
(1993). 

+oo 30. L. Onsager, Nuovo Cimento Suppl. 6, 279 (1949). 

G _ y _ T" 2 / In (T + e~ x ) dx P-H- Chavanis and J. Sommeria, J. Fluid Mech 356, 259 

Jo K ' (1998). 

H 82) 32 - B - Turkington, A. Majda, K. Haven, M. Dibattista, PNAS 
98, 12346 (2001). 

U sing 33. P.H. Chavanis, Physica D 200, 257 (2005). 

34. R. Robert and J. Sommeria, Phys. Rev. Lett. 69, 2776 

r+°o n 2 (1992). 

/ In (l + e~ x ) dx = — , (183) 35. P.H. Chavanis, Phys. Rev. E 65, 056123 (2002). 
Jo ±* 

we establish Eq. ijS^I. 
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